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The  overestiaation  of  power  (tensity  by  ths  Fraunhofer 
( "far-field" )  approximation  to  ths  more  exact  Fresnel  ( "near-field" ) 
diffraction  integral  is  determined  as  a  function  of  Fresnel  muster 
at  arbitrary  on-axis  field  points  for  rectangular  and  circular 
aperture  illuminations  of  unifora  amplitude  and  phase.  For  square 
and  circular  apertures,  the  overestiaation  is  2.8  percent  and  1.3 
percent,  respectively,  at  a  Fresnel  number  of  1/8  ( corresponding  to 
the  far-field  boundary  distance  2D*/M  and  518  percent  and  147 
percent,  respectively,  at  a  Fresnel  mater  of  1. 


♦This  work  was  sponsored  by  the  Electronic  Systems  Division,  Air 
Force  Systeas  Command,  Hanscom  Air  Force  Base,  Massachusetts. 
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SECTION  1 


INTRODUCTION 


The  Fraunhofer  and  Fresnel  diffraction  integrals  are  coononly 
used  analytical  tools  for  the  evaluation  of  fields  originating  from 
illuminated  apertures.  These  integrals  are  derived  in  appendix  A 
from  the  Fresnel-Kirchhoff  integral  subject  to  small  angle 
conditions.  Whereas  the  exponential  factor  of  the  Fresnel  integrand 
retains  both  linear  and  quadratic  phase  terms  (as  a  function  of  the 
aperture  coordinates)  and  is  applicable  to  near-field  points  as  well 
as  far-field  points,  the  Fraunhofer  integrand  retains  only  linear 
phase  terms  and  is  applicable  only  to  far-field  points.  The 
Fraunhofer  "far-field"  formulation  is,  therefore,  an  approximation 
to  the  more  exact  Fresnel  "near-field"  formulation.  (This  statement 
is  not  applicable  at  a  lens  focal  point  where  both  linear  and 
quadrature  phase  terms  are  zero.) 

For  on-axis  field  points  and  an  aperture  illumination  of 
uniform  phase,  the  Fraunhofer  phase  terms  are  zero,  corresponding  to 
totally  constructive  interference.  However,  for  the  same 
conditions,  the  Fresnel  phase  terms  are  non-zero,  corresponding  to 
destructive  interference  among  interfering  rays.  Consequently,  the 
far-field  formulation  yields  a  larger  power  density  at  an  arbitrary 
on-axis  field  point  than  the  more  exact  near-field  formulation  for 
an  aperture  illumination  of  uniform  phase. 

In  this  paper,  the  overestimation  of  power  density  by  the 
far-field  formulation  is  determined  as  a  function  of  Fresnel  number 
at  arbitrary  on-axis  field  points  for  rectangular  and  circular 
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aperture  illuminations  of  uniform  amplitude  and  phase.  For  square 

and  circular  apertures ,  the  overestimation  is  found  to  be 

2.8  percent  and  1.3  percent,  respectively,  at  a  Fresnel  number  of 

2 

1/8  (corresponding  to  the  far-field  boundary  distance  2D  /X)  and  518 
percent  and  147  percent,  respectively,  at  a  Fresnel  number  of  1. 

The  1.3  percent  result  is  compatible  with  a  result  reported 


earlier 


(1) 


The  Fresnel  number  is  defined  in  section  II.  Far-field  and 
near-field  formulations  are  summarized  in  sections  III  and  IV, 
respectively.  A  comparison  of  far-field  and  near-field  formulations 
is  given  in  section  V. 


SECTION  2 


DEFINITION  OF  FRESNEL  NUMBER 


Consider  a  collimated  beam,  of  uniform  amplitude  and  phase  with 

2 

intensity  I0  (W/m  )  at  a  wavelength  X(m) ,  which  is  incident  on  a 

limiting  aperture.  For  a  rectangular  aperture,  the  clear  aperture 

is  of  width  2b  and  height  2b  (see  figure  la).  For  a  circular 
x  y 

aperture,  the  clear  aperture  is  of  radius  b  (see  figure  lb).  In 
figure  1,  a  collimated  beam  corresponds  to  a  source  point  S  at  a 
distance  d'  -  •. 

For  a  circular  aperture  of  radius  b,  the  number  N  of  Fresnel 
zones,  subtended  by  every  point  on  the  edge  of  the  circular  aperture 
at  a  field  point  P(0,0,d)  on  the  optical  axis  at  a  distance  d  from 
the  aperture,  is  defined  as 

N  ■  l***')^**6')  «  b2/XL,«  b«d,  b«d'  (1) 

where 

s  -  (b2  +  d2)1</2,  s'  *  (b2  +  d'2)1^2  (see  figure  lb) 

L  ■  [(1/d)  ♦  (1/d' ) J”1 

For  a  collimated  beam,  L  -  d  and  equation  (1)  reduces  to 

N  «  b2/Xd;  b«d,  collimated  beam  (2) 

For  points  P  in  the  very  near  field,  N  ■*  •.  For  points  P  in 

the  very  far  field,  N  ■+  0.  For  a  collimated  beam,  the  boundary 

distance  between  the  near  field  and  far  field  is  conventionally 
(1)  2 

chosen  asv  d  -  2(2b)  /X  which  corresponds  to  a  Fresnel  number 
N  -  b2/Xd  -  1/8. 


ting  Aperture 
Rectangular  A 
Circular  Aper 


Consider  now  the  rectangular  aperture  of  figure  la.  Not  every 
point  on  the  edge  of  the  rectangular  aperture  subtends  the  sane 
number  of  Fresnel  zones  (except  in  the  very  far  field). 

Consequently,  the  effect  of  edge  diffraction,  on  the  intensity  at  a 
near-field  point,  is  less  pronounced  for  a  rectangular  or  square 
aperture  than  for  a  circular  aperture.  Edge  diffraction  by  a 
rectangular  aperture  is  characterized  by  the  Fresnel  numbers,  Nx  and 
Ny,  defined  by 

sx  -  d  2 

Nx  ”  “v2 —  58  bx  bx«d,  collimat®d  beam  (3a) 

sv  -  d  2 

Ny  ”  X/2 —  “  y  '  by«d,  collimated  beam  (3b) 

where 

sx  "  (bx  +  d2)1/2'  sy  *  (by  +  d2)1/2  (see  figure  la) 

The  Fresnel  number  Nx  is  the  number  of  Fresnel  zones  subtended  by 

every  point  on  the  edges  x  -  +  bx  of  the  rectangular  aperture  at  the 

field  point  P  if  the  point  source  S  were  replaced  by  a  line  source 

of  infinite  extent  parallel  to  the  y  axis  and  if  the  aperture  were 

of  infinite  extent  in  the  y  direction.  A  similar  statement  applies 

to  N  if  x  and  y  are  interchanged  in  the  statement  for  N I. 
y 


SECTION  3 


FAR- FI ELD  FORMULATION 


2 

The  on-axis  intensity  I  (W/m)  at  a  far-field  point 

P(x,y,z)  -  P(0,0,d)  is  found  from  the  Fraunhofer  diffraction 

12) 

integral  of  appendix  A  to  be'  ' 


1(0, 0,d)  -  |u(0,0,d)|2  -  J  J  dxDdy( 


xoA 
A2  d2 


PtA 
A2  d2 


Pt* 

4nd2 


,  collimated  beam 


(4) 


where 


A  -  aperture  area 


{ 


4b  b  r  rectangular  aperture 
X  y 


lib 


circular  aperture 


Pt  -  total  power  radiated  -  IQ  A 

4njk 

g  -  aperture  (antenna)  gain  ■  * 

( X/2b  HV&  )  *  rectangular  aperture 


U  n/(  A/2b)  ]' 


,  circular  aperture 


It  will  be  noted  from  equation  (4)  that  the  beam  is  no  longer 
collimated  in  the  far  field  but  appears  to  be  spreading  as  though 
the  aperture  were  a  point  source  whose  radiation  is  confined 
primarily  to  beam  total  angles  of  the  order  of  (A/2bx)  in  the  x-z 
plane  and  (A/2by)  in  the  y-z  plane. 
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The  intensity  I  at  a  field  point  P(0,0,d)  can  be  expressed  in 
terns  of  the  number  of  Fresnel  zones  subtended  by  the  aperture  at 
the  field  point.  Noting  from  equation  (2)  and  equation  (3)  that 
A  -  ixNXd  for  a  circular  aperture  and  that  A  -  4-Jn^  'ITT  Xd  for  a 
rectangular  aperture,  equation  (4)  reduces  to 

'  lQ  16  Nx  Ny  ,  rectangular  aperture  (5a) 

1  “  2  _2  ' 

I  »r  fr  ,  circular  aperture  (5b) 

o 

collimated  beam,  N  ,  N  ,  N  <  1/8 

x  y 

It  will  be  noted  from  equation  (4)  or  equation  (5)  that  in  the 

2 

far-field  I  -*  0  as  1/d  -»  0  or  equivalently  I  ■+  0  as 

N  N  (or  N2)  -»  0. 
a  y 


SECTION  4 


NBMt-PXELD  rOMULKTIOM 

The  intensity  I  (W/k2)  at  any  on-axis  fisld  point  p(0,0,d)  of  a 
rectangular  aperture  with  illuadnation  of  uni  fora  amplitude  and 
phase,  is  found  fran  the  Presnel  integral  to  be*3* 

I  -  4IC  [(^(/S^)  ♦  [c2(^T^)  ♦  «2(/3^)]»  (6) 

coll  ina ted  bean,  bx«d,  by«d 

where 

w 

C(w)  -  J  cos[(x/2)t2]  dt  •  Preenel  cosine  integral * 4 * 

0 

w 

S(w)  -  J  sin£(s/2)t2j  dt  -  Preenel  sine  integral*4* 

0 

Nx,  Ny  -  Presnel  nusfcers  defined  by  equation  (3) 
lQ  -  aperture  intensity  (H/la2) 

The  Presnel  integrals  have  the  properties  that*7* 


S(w)  - 

-  S(-w),  C(w)  -  -  C(-w) 

(7a) 

SC)  - 

CC)  -  % 

(7b) 

S(0)  - 

C(0)  -  0 

(7c) 
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rot  Nx  -  Ny  -  equation  (6)  reduces  to  I  ■  IQ  which  is  to  be 
expected  for  points  in  the  very  near  field.  For  Mm  •  •  0, 

equation  (6)  reduces  to  I  -  0  which  is  to  be  expected  for  points  in 
the  very  far  field. 

Equation  (6)  reduces  to  equation  (Sa)  in  the  Unit  of  Mall 

/C) 

Fresnel  nuabers.  this  result  is  rtaarmst  rated  by  noting  that1 


C(w) 

-  w  - 

♦  ••• 

(6a) 

S(w) 

.3  .  .... 

(6b) 

c2 

(w)  ♦  S2(w) 

-w2- 

'  TS  ♦  •••• 

(6c) 

♦  32(/arp 

■*. 

-  (*/2>2(av3  ♦  .... 

(9a) 

c^/asj) 

♦  s2</a^) 

-** 

-  (V2)2(2hy)3  ♦  .... 

(9b) 

Substituting  equations  (9a)  and  (9b)  with  the  condition  (4e2/45)(M^ 

♦  nJ)  «  1,  into  equation  (6),  I  •  IQ  16  which  is  the  result 

given  by  equation  S(a). 

For  a  circular  aperture  which  subtends  a  snail  angle  at  a  field 
point  F(0,0,d) ,  the  on-axis  intensity  X  is  found  froa  the  Fresnel 
integral  to  be 


I 


~  J  «p  t  dr  p 

0 

-  I0  I  1  -  np  (-!•#)  I2-  Ie  {(1  -  cos(*))2*  sin2( XN)} 

•  4Iq  sin2(«W/2);  col  11m tad  bean,  b«d  (10) 

where 

r  -  (x2  ♦  y2)172 
N  -  b2/Xd. 


ror  N  -  1,  3,  5,  ....,  I  -  4Iq.  Por  W  -  2,  4,  6,  I  -  0. 

for  n  -  I  •  zQ.  (aquation  (10)  does  not  extrapolate  to  this 
rasult  because  aquation  (10)  ia  not  valid  Cor  apertures  which 
subtend  a  large  angle  at  the  field  point.  See  reference  (6)  for 
circular  apertures  which  subtend  a  large  angle  at  the  field  point.] 
with  the  substitution  exp  (-i«)  *  1-iMN  for  (*V/12)  «  1 , 
aviation  (10)  reduces  to  I  *  l0x2w2  which  is  coapatible  with  the 
tar-field  result  given  by  equation  5(b). 
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SECTION  5 


COMPARISON  OP  PAR-FIELD  AND  NEAR-FIELD  POMULATIONS 


The  on-axis  near-field  intensities  for  square  apertures  of  side 

length  2b  -  2b  -  2b  and  circular  apertures  of  disaster  D  -  2b  with 

*  y 

illumination  of  uniform  amplitude  and  phase,  are  compered  in 
figure  2  as  a  function  of  the  Fresnel  number  N^  -  Ny  ■  N  -  b2/Xd. 

The  near-field  formulations  given  by  equations  (6)  and  (10)  are  also 
compared  with  those  obtained  from  the  far-field  formulations  given 
by  equations  (5a)  and  (5b). 

The  peak  intensities  attained  are  4IQ  at  Fresnel  nusbers  M-l, 

3,  5,  ...  for  circular  apertures  and  3.24  IQ  at  a  Fresnel  number 
N  •  3/4  for  square  apertures.  The  minimum  intensities  attained  in 
the  near  field  are  0  IQ  at  Fresnel  nusbers  N  •  2,  4,  0  ...  for 
circular  apertures  and  0.37  lQ  at  a  Fresnel  nuaber  N  -  7/4  for 
square  apertures.  The  intensities  obtained  from  the  far-field 
formulations  exceed  those  obtained  from  the  more  exact  near-field 
formulations. 

The  ratios  of  these  intensities  are  given  in  table  1  as  a 
function  of  the  Fresnel  nuaber  N.  For  square  and  circular 
apertures,  the  far-field  formulation  overestimates  the  power  dsns 
by  2.8  percent  and  1.3  percent,  respectively,  at  a  Fresnel  nuaber 
1/8  (corresponding  to  the  far-field  boundary  distanoe  ID2/!)  and 
percent  and  147  percent,  respectively,  at  a  Fresnel  number  of  1. 

The  overestimation  (d84/2)2/sln2(i84/2))  -  1  for  circular  apertures 
has  bean  reportsd  earlier^ * .  Please  note  the  round-off  error  in 
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reference  (1)  for  a  Fresnel  nuaber  of  1/4  (G/GQ  a  0.95  and  not 
0.94).  For  Fresnel  nuabers  corresponding  to  constructive  or  weakly 
destructive  interference,  the  overestiaation  is  greater  with  square 
apertures  than  with  circular  apertures.  However,  for  Fresnel 
numbers  corresponding  to  totally  or  almost  totally  destructive 
interference,  the  overestiaation  is  greater  with  circular  apertures 
(see  N-2  and  N-2.25  in  table  1).  Constructive  and  destructive 
interference  are  more  pronounced  for  circular  apertures  and, 
therefore,  are  in  better  and  worse  agreement,  respectively,  with  the 
totally  constructive  on- ax is  interference  of  the  Fraunhofer 
integral. 
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appendix  a 


DERIVATION  OF  THE  FRAUNHOFER  AND  FRESNEL  DIFFRACTION  INTEGRALS 

For  a  source  point  S(x',y',d'),  arbitrary  aperture  point 
q(x  ,yG)  with  origin  0,  and  field  point  P(x,y,d),  the 
Fresnel-Ki rchhof f  diffraction  integral  for  the  amplitude  U(P) 
(proportional  to  either  the  scalar  electric  or  magnetic  field 
intensities)  is  given  in  the  notation  of  figure  1  by 

U(P,  -  -  fi  JJ  Ico»(S,S')  -  cos  <i«»d*0dy0  <*-!> 

A 

where 

s'  -  §3,  s'  -  | SQ | 
t  -  p3,  s  -  |P$| 

n  -  unit  vector  normal  to  the  aperture  in  the  direction  of  the  +z 
axis 

A  -  area  of  the  clear  aperture 

B  -  complex  constant  ■  d'  /T~  e  ^  which  satisfies  the  bou  dary 
condition  U(0,0,0)  »  . 

(The  aperture  point  Q(xq,  yQ)  is  not  restricted  to  the  aperture 
point  Q  shown  in  figure  1). 

For  the  small  angle  conditions 

(n,s' )  »  (n,  §3),  (n,s)  «  (n,P<5)  (A^2) 

and  the  straight  line  condition 
(§5,  p3)  »  r  radians 

(1)  H.  Born  and  E.  Wolf,  "Principles  of  Optics  "(Pergamon  Press, 
Oxford  1964),  pp.  382-383. 


A- 1 


then 


cos(n,s')  a  1,  cos(n,s)  a  -l,  [conditions  (A-2)  and  (A— 3) )  (A-4) 


1/s's  a  1/d'd  [condition  (A-2)] 

Expanding  s'  and  s  in  a  power  series, 

-  ,  .  ,  (xxo  *  Wo>  ,  *0  *  4  <xxo  *  Wo>2 

8  +s’d  td - a - +  — - ^ - • 

exp  (ik(s'+s)]  -  exp  [ik(d'+d)J  exp  (ikf) 
where 

(*y  vy0)  .  i  xp  +  4  <**<,  *  w0»2 

d  5  I.  '2d3 

L  «  [(1/d)  +  (l/d')]_1 

Substituting  equations  (A-4)  -  (A-7)  and  the  boundary  condition 
U( 0,0,0)  -  /F^  into  equation  (A-ll), 


(A-5) 

•  (A-6) 
(A-7) 


~  i/rn  rr 

U(P)  -  — 5-3 —  exp  (ikd)  JJ  exp(ikf)  dx  dy 
au  A 


o  ■‘o 


(A-8) 


If  only  the  first  term  in  the  above  expansion  for  f  is  retained, 
then  equation  (A-8)  reduces  to 


-  i/T 


U(P)  -  ^  2  exp  (ikd)  JJ  exp  [(-ik/d)  (xxQ+  yyQ)]  dxQdy0  (A-9) 


Equation  (A-9)  is  known  as  the  Fraunhofer  diffraction  integral . 

If  only  the  first  two  terms  in  the  above  expansion  for  f  are 
retained,  then  equation  (A-8)  reduces  to 

U(P)  -  ^  ~  exp  (ikd)  (JJ  exp  [(-ik/d)(xxQ+  yy0)] 

A 

exp  [(ik/2L)(x*+  y*)]  dxQdy0  )  (A-10) 

Equation  (A-10)  is  known  as  the  Fresnel  diffraction  integral. 


A- 3 


p  jtn  e 


n?7 


•  -  • 


